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Thermal fluctuations are important for amphiphilic bilayer membranes since typical 
bending stiffnesses can be a few k B T. The rod-like constituent molecules are generically 
tilted with respect to the local normal for packing reasons. We study the effects of fluc- 
tuations on membranes with nematic order, a simplified idealization with the same novel 
features as realistic tilt order. We find that nematic membranes lie in the same universality 
class as hexatic membranes, i.e. the couplings that distinguish nematic from hexatic or- 
der are marginally irrelevant. Our calculation also illustrates the advantages of conformal 
gauge, which brings great conceptual and technical simplifications compared to the more 
popular Monge gauge. 
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The study of fluctuating surfaces in an intriguing statistical mechanics problem with 
numerous applications. Amphiphilic molecules in water can self-assemble into thin flexible 
bilayer membranes which provide an experimental realization of random surfaces . As we 
recall below, a membrane's internal order determines the nature of its shape fluctuations. 
Although fluctuating membranes with hexatic bond-orientational order have received a 
lot of theoretical attention @ [[§ [H , the generic situation is for the rod-like constituent 
molecules to tilt with respect to the local surface normal. We shall see that tilt order 
differs from hexatic order since tilt allows certain anisotropic couplings. In this letter we 
consider "nematic" membranes, the simplest membrane model with in-plane orientational 
order and anisotropy. We find the anisotropy to be marginally irrelevant; these membranes 
may be considered hexatic at very long length scales. 

Our framework is continuum elastic theory since the micron size of membrane struc- 
tures is much larger than the typical molecular size, a few tens of Angstroms. If the 
molecules are allowed to take up their preferred area, then surface tension will be unim- 
portant and bending rigidity will dominate. Typical bending moduli for lipid bilayer 
membranes are on the order of 10k B T ||, which is low enough for thermal fluctuations to 
play a role. We study fluid-like membranes in which the molecules are free to diffuse along 
the membrane in response to shape fluctuations. Therefore shape and orientational order 
are the only elastic degrees of freedom, and our elastic theory must be coordinate-invariant 
i- 

The simplest type of in-plane orientational order is hexatic order H- Pure hexatic 



order is readily seen in Langmuir monolayers |§ and thin smectic liquid crystal films ||10|| , 
but this order has not to our knowledge been confirmed in bilayer membranes. Instead, 
in the ordered phases the constituent molecules generically tilt for packing reasons. If tilt 



and hexatic order are present together [pT| , then in the absence of defects we can for the 
purposes of elastic theory treat these as locked together. Tilt and hexatic order break the 
same rotational symmetry, and together lead to just one elastic mode. Since tilt order has 
less symmetry than hexatic order (see below), the relevant order is tilt when both tilt and 
hexatic order are present. 

One can represent hexatic order on a membrane by a unit tangent vector field 0. 
Each vector points from a molecule to one of the molecule's nearest neighbors and is 
defined up to 2tc/6 rotations. Demanding this symmetry automatically brings along a 
larger symmetry, i.e. global rotations through an arbitrary angle; in this sense hexatic 
order is isotropic. This is reminiscent of the isotropy of the continuum elastic theory of a 
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two-dimensional triangular lattice [T!|, but holds even on curved surfaces, whose principal 
curvatures could in principle have cared about a global rotation of in-plane order. 

Similarly one can represent tilt order as a unit tangent vector field rh defined by the 
direction of the projection of the axis of each molecule on the local tangent plane. We 
consider only the elastic modes so we fix the polar angle between the normal n and the 
molecular axis. Bilayers lack a preferred normal, so to complete our specification of rh 
we demand the bilayer symmetry n — > —h, rh — > —rh ||13|| ||14| . In contrast to the hexatic 



case, the energy for bending a membrane along an axis parallel to the direction of tilt can 
clearly be different from the energy to bend the membrane the other way. (Mathemati- 
cally we will see in a moment how this conclusion arises from the reduced symmetry of tilt 
order.) This tilt-shape coupling has a number of simple consequences; for example, tilt 
order is responsible for the Pat ripple phases [|TT||15|]. Tilt order is also manifested dra- 
matically when the constituent molecules are chiral: the microscopic molecular chirality 



affects macroscopic membrane shape only in the presence of tilt [16|. Tilt order can play 
a role in the budding of artificially manufactured vesicles [|T7] and it may be the typical 
order in biological membranes [TOl . 

The nature of membrane shape fluctuations depends on the degree of internal order. 
We characterize these fluctuations by using the Wilson renormalization group to compute 
the long-distance behavior of the elastic couplings. For example, the bending rigidity of a 
fluid membrane becomes ineffective beyond a persistence length £ P , leading to a crumpled 
phase at any nonzero temperature [|n|.0 Internal order tends to stiffen a membrane. For 
stiff enough elastic constants, the bending rigidity for a hexatic membrane reaches a fixed 
point at very long length scales, possibly leading to a "crinkled" phase with quasi-long- 
range order in the normals 0. Finally, non-self avoiding tethered membranes (crystalline 
membranes with an infinite core energy for dislocations) undergo a crumpling transition 
and have a flat phase at low temperature [ffijf . It is natural to extend these analyses to 
the case of tilt order, since as we remarked it is experimentally the most relevant regime 
with orientational order. To keep our formulas compact we will impose a "nematic," or 
2-atic symmetry to get the simplest model with anisotropy; the same technique applies to 
the more realistic case with no discrete symmetry. We will then ask, does the anisotropy 
lead to new fixed points and new physics, or is the anisotropy marginally irrelevant as in 
flat thin liquid crystal films H? 



1 We ignore self-avoidance since we will ultimately study a fixed point at large stiffness. 
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To answer these questions with precise calculations, we must choose coordinates. 
The most popular choice for such calculations is Monge gauge, in which the surface is 



parametrized by its height above some flat reference surface |T!|[]14]]. It turns out that 



Monge gauge is not very convenient for tilted membranes; conformal gauge, as used for 



example by Polyakov ||21]1 , is much better-suited. Conformal gauge has the technical ad- 
vantages of leading to compact expressions and thus less algebra than Monge gauge, as 
well as the conceptual advantage of making explicit a useful separation between extrinsic 
and intrinsic geometry. 

We begin by describing the geometrical constructions necessary to write down the tilt 
free energy. Even with the benefits of conformal gauge, the formulas can get a bit long. 
Thus, we make some inessential but simplifying assumptions. Next, we briefly summarize 
the issues involved in conformal gauge calculations. Finally, we present our result and 
discuss its implications. 

Following the notation of ||21|| . we parametrize our surface by x(£), where if is a three 
dimensional vector and (^ 1 ,^ 2 ) are the (arbitrary) two dimensional coordinates. From 
x(£) we construct the metric tensor g ab = d a x ■ d b x and the second fundamental form 
K ab = n ■ d a d b x, where d a = d/d^ a and n is the local normal. We denote the covariant 
derivative associated with g ab by V, the inverse of g ab by g ab , and the determinant of g ab 
by g. It is useful to introduce a local set of orthonormal frames e a ; we denote by e a a 
the change-of-basis matrix that converts orthonormal frame indices ... to coordinate 
indices a, 6, .... Thus for example V a m a = <9 a m a +e a /3 fl a m^ where fi a is the spin connection 

I- 

The elastic free energy for tilted bilayer membranes must have coordinate invariance, 
Euclidean invariance, and the discrete bilayer symmetry discussed above. Writing only the 
bulk terms quadratic in the curvature or derivatives of m, i.e. only the marginal terms, 
we find F tilt = F 1 + F 2 where 

Fx =\ j d 2 ^[^i(V • rh) 2 + k 2 (V x m) 2 + K(K a a ) 2 

+ K x rh ■ K ■ K ■ rh + K 2 (m ■ K ■ rh) 2 ] 
F 2 =\ j d 2 ^[^-K-m)(V-m)+(3 2 (K a a )(V-m) 

+ P 3 K\V a m b + [3 4 m c K ca m b V b m a ] (1). 



In ([!]), we contract indices with g ab , and V x rh = e a b V a m b , where e ab is the covariant anti- 
symmetric tensor .i Note that F x reduces to the hexatic membrane free energy when k 1 = 
k 2 = and k 1 = k 2 = k A , since even on a curved surface (V-m) 2 + (V x m) 2 = V a m b 'V a m b . 
Following Nelson and Pelcovits || , we find it convenient to take k 2 > k 1 without loss of gen- 
erality and rewrite the terms involving only m and the metric as k 2 V a m b V a m h + k{V ■ rh) 2 , 
where k = k 2 — k 1 . 

As we alluded to above, for large stiffnesses k, k A , hexatic membranes are governed by 
a line of fixed points in the K~ x —k A plane [[§. Our first simplification for tilt is to study 
the stability of this line against anisotropy by working to first order in the anisotropic 
couplings. We further simplify our job by dropping the terms of F 2 in ([!]). This truncation 
is mathematically consistent because F 1 is a complete list of terms with the "nematic" 
rh — > — rh symmetry. We do not know of an experimental system with this symmetry. One 
can imagine a membrane made of rod-like molecules that lie parallel to the local tangent 
plane, or perhaps more realistically a membrane with stiff rod-like molecules aligned along 
the local normal but with a rectangular cross-section. In any case, this free energy is 
the simplest membrane model with in-plane order and anisotropy; we expect it to be 
qualitatively similar to the case of real tilt. 

We study the long-distance behavior of our model with the Wilson momentum-shell 
renormalization group |2^] . The first step of this procedure is to decompose the fields into 



slowly varying and rapidly varying parts x = x + x 1 ,rh = rh Q + fh 1 , where x ,rh have 
Fourier modes with wavevector k satisfying \k\ < A/6, and x x , fh 1 have Fourier modes in a 
shell A/6 < <A. A = 27r/ais the wavevector cutoff corresponding to the short-distance 
cutoff a, and b is a number slightly greater than one. The next step is to coarse-grain the 
system by tracing over the fast modes in the partition function. We write the partition 
function for the membrane as a path integral over the shape and orientational degrees of 
freedom Since membranes are stiff but not completely rigid, we work to first order 
in T/K.i This approximation amounts to integrating out the fast modes in the Gaussian 
approximation. Finally, we rescale the distances to restore the cutoff to its original value 
and obtain the recursion relations for the effective couplings. 



2 In (|l|) we have corrected a redundancy in [14]; we note here that the a.\ and a.2 terms of [14] 
differ by a total derivative. 

3 Since we work in the stiff regime, we ignore defects in the in-plane order. These will not 
qualitatively change our result (but see [S]). 



To carry out these calculations, we must choose coordinates. The metric g ab = d a x-d b x 
at each point of a two dimensional surface is a symmetric 2x2 matrix and thus has 
three independent degrees of freedom. Two of these are removed by the two independent 
coordinate degrees of freedom, leaving one physical degree of freedom. In fact, we can 
always find local coordinates in which the metric is a spatially dependent conformal factor 
times the trivial metric M: 

9 ab (0 = P(0<U- (2) 

Conformal gauge is the two dimensional analog of arc-length parameterization of a curve 
in space. 

The tensor fields associated with a surface have a compact form in conformal gauge. 
For example, the mean curvature is ft ■ V 2 x = p~ x h ■ d 2 x, the change of basis matrix is 
e a a = VP~ S aa, and the s P in connection is fi a = ±e ab d b \ogp 0. 

For purposes of illustration we first consider the fluid membrane (no in-plane order) in 
conformal gauge, as discussed by Polyakov |23J] . Since p depends on x, the fluid membrane 
free energy is a complicated nonlinear functional of x and not in the most convenient form 
to carry out the renormalization group procedure. To treat p as an independent field, we 
introduce a delta function constraint and then enforce the constraint by introducing some 
Lagrange multiplier fields [pl |: 



[dx]exp [-A J d^p'^x) 
[dx\ [dp}5[d a x ■ d b x - p5 ab ] exp 
[df][dp][dA a6 ]e- F 



1/q2 -?\2 



where 



F= 2^/ ^^P~ l ^ 2 x) 2 + t\ a \d a x-d h x-p5 ah 



(3) 



(4) 



To complete the specification of the partition function, we must define the functional 
measure. There are various factors that arise from coordinate invariance [^3[. We simply 
note that the geometrical measure factors will not enter our calculations to 0(T/k), and 
the Liouville counterterm will not enter the recursion relations for the bending stiffnesses. 
A further advantage of conformal gauge over Monge gauge is that these conformal gauge 
measure factors have been well-studied [^3| , making conformal gauge the better choice for 
calculations beyond 0(T/k). 
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The Wilson renormalization group procedure requires the most general low order ex- 
pression consistent with all the symmetries; ([|) is not obviously of this form. The Lagrange 
multiplier field enters in a very specific way. More precisely, since the renormalization group 
is an iterative procedure, we must be sure that after eliminating short wavelength modes 
the long-wavelength effective free energy differs from the original only by the values of the 
couplings. If e.g. a A 2 term were generated, then the delta function constraint would be 
softened to a Gaussian and we could not define the recursion relations. Such terms do in 
fact appear, but always suppressed by powers of the short-distance cutoff. 

The extra fields X ab and p were supposedly introduced to simplify the calculation; 
it is natural to ask why the extra fields do not lead to more complexity. The continued 
sharpness of the delta function during coarse-graining leads to a technical simplification. 
To see this, we recall Polyakov's trick to diagonalize the fluctuation part of the free energy 



functional |2l]]. He decomposes the rapidly varying part A" of the Lagrange multiplier 



field X ab into a traceless and transverse trace part taking 

Af = dj lb + d b f la - S ab dJ lc + (d ab - Ci- (5) 

Using this decomposition, we expand the free energy to quadratic order in the fast fields 
and split its quadratic part F H into a free part and a part to be treated as a perturbation 



^ = F A +F, 



B ! 



f a = ^J d^fof 1 ^) 2 - 2t(d 2 f la )(d a x ) ■ x, - i( lPl ] 



(6) 



F B ~ J d^atfidj,) ■ {d b x x ) - 2p lPo - 2 (d 2 x ) ■ (d%) 
+ p 1 2 p Q - s (d 2 x ) 2 -2 i \ a 1 \d a d b x )-x 1 ]. 

In the diagrammatic expansion of the effective free energy, F A determines the propagators 
and F B determines the vertices. Since there is no d 2 term in (|||), the p x propagator is 
zero: (piPi) = 0. Many diagrams that otherwise would have contributed are thus zero. 
Also, the bending free energy for a fluid membrane is independent of derivatives of p, so 
we can treat p Q as a constant. Thus the conformal factor p does not play much of a role 
in the calculations. 

Turning to our model, we introduce another Lagrange multiplier p to enforce rh 2 = 1. 
Introducing p is more convenient than parametrizing rh as m a = (cos 9, sin 9) a , since all the 
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terms of ([!]) are nonlinear in 9. Working directly in terms of fa is covariant and makes for 
compact expressions. Also, the Lagrange multiplier trick takes care of field renormalization 
automatically (cf. [0). 

We are now ready to outline our calculation. Denoting the fields (x, p, /, £, fh, /x) 
collectively as (f), we want to compute the effective free energy 



log 



/ [#i] 



exp 



1 



TrlogO[4> } 



(8) 



where O[<p ] is the matrix associated with the fluctuation free energy F n = | j d 2 ^(j) 1 O[(p ](p 1 . 
We split O[4> ] into a free part A [</> O ] and an interaction part O B [4> }. Since we work to 
first order in the anisotropic couplings fc, k 1 , k 2 , we treat all the subterms in the expansion 
of the anisotropic terms as perturbations. Thus, 



O 



A 



( Po^Xb 

-id 2 (d a x ) 



V 



-i(d b x )d 2 




\ 



2fl2 



(9) 



/ 



Here we can see that the block diagonal form of O a (and of O a _1 , the matrix of 
propagators) leads to a clear separation between extrinsic and intrinsic geometry. There 
are no propagators connecting x x fields (extrinsic) with p x or rh 1 fields (intrinsic). In our 
expansion, the only coupling between extrinsic and intrinsic geometry is through graph 
vertices derived from the Lagrange multiplier term of (|j). 

The free energy is given by expanding the logarithm in (§) in powers of O a ~ 1 O b : 



\^\og{0 A + O b ) = \Tr[logO A + O a 1 O b - \o a 1 O b O a 1 O b + 



(10) 



The first term of ( [10]) corresponds to graphs with no vertices, the second to graphs with 
one vertex, and so on. Since O a depends on the slow fields we must keep the first term; 
unlike in Monge gauge, graphs with no vertices contribute to the renormalization of the 
stiffnesses. Evaluating the traces of ([TOD leads to the effective free energy 



eff 



1 

2T 



2tvk 



3 5 £ 17 
~ K+ 8 2+ 16 + T Kl 
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- K, 



log b 



+ 



2T 



T 



k± H (—k — 7k 1 + 11k 2 ) l°gk 



d 2 ^Po lfn o (d a d a x ) ■ (d a d x o )m% 



+ 



1 



1 

+ 2T 



2T 



K 2 



Tk 2 
8tyk 



log b 



d 2 £p lm O m O m O m o( d a d f3^o) ■ ( d -y d 5^o) 



Po 



ab 



~8k 



m 0a m 0b 



^1 — ^" ^2 

2/T~ 



))log& 



+ i / d-«^)> (i + ^icg.) + A / d^(«) 2 (i - j-io g *) . (ii) 



In deriving (P|) , we have used 9 a x • 9 b x = p 5 ab + 0(T/k) and 



l + 0(T/«). 



The purpose of quoting the long expression (|TT|) is to point out an unexpected term: 
the term with X ab m 0a m ob did not appear in the original free energy. This term spoils the 
conformal gauge condition d a x ■ d b x Q = p S ab for the slow fields, but it does not spoil the 
renormalization group calculation. We simply integrate over Aq & to get a modified delta 
function relating d a x -d b x , p S ab , and m 0a m ob . Using the delta function to integrate over 
p , we obtain a free energy of the same form as (JJ) but with modified coefficients. Sub- 
stituting the hexatic fixed line relation k 2 = 4k, we find the linearized recursion relations 
dg/dlog& = Mg, with 
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/ Bk 8k 



V 



(12) 



g = ((/3k) -1 , (/3fc 2 ) _1 , /3A;, /?K l5 /3k 2 ), an d = V^- The matrix M has eigenvalues 
— 3/(47T/9 2 k 2 ) and corresponding to the hexatic membrane, and three more negative 
eigenvalues, (— 9 ± -v/41)/(16tt/3k), — 3/(8tt/3k), corresponding to the anisotropic couplings. 

The hexatic fixed line is therefore stable to the anisotropic terms obeying the rh — > —rh 
symmetry. Beyond a length scale exponential in the stiffness k/T we expect nematic 
membranes to behave as hexatic membranes. This is in accord with our intuition that 
since there is no true long-range in-plane order in this two dimensional system, at very 
long length scales anisotropy should be irrelevant. The exponential dependence of these 
length scales on the stiffness means in practice anisotropic membranes will differ from their 
hexatic counterparts. We remark that Peliti and Prost have also argued that anisotropic 
membranes should lie in the same universality class as hexatic membranes |T3[]. However 
our calculation shows that the actual way in which this happens is different from the 
screening mechanism they proposed. 

We have seen that conformal gauge is well-suited for determining the long-distance 
properties of fluctuating membranes with in-plane order. We can easily include the terms 
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we ignored to extend our results to more realistic tilted membranes. Our technique is 
useful for other applications. One example is the critical exponents of a crinkled tilted 
membrane; we expect these exponents to differ from the hexatic case since the anisotropy 
is only marginally irrelevant (c/. ||). 
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